Water wave collapses over quasi-one-dimensional non-uniformly periodic bed profiles 
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Nonlinear water waves interacting with quasi-one-dimensional, non-uniformly periodic bed pro- 
files are studied numerically in the deep-water regime with the help of approximate equations for 
envelopes of the forward and backward waves. Spontaneous formation of localized two-dimensional 
wave structures is observed in the numerical experiments, which looks essentially as a wave collapse. 
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Nonlinear waves of different nature, when propagating 
in spatially periodic media, are known to introduce var- 
ious interesting phenomena. In particular, the so called 
gap soiitons (GSs) can be mentioned, which are self- 
localized waves existing due to the presence of a Bragg 
gap in the spectrum of linear waves, and, on the other 
hand, due to nonlinear interactions between components 
of the wave field. Gap soiitons are studied mainly in the 
nonlinear optics (see, e.g., Refs.pl H H II H HE H !, 
llOlllllrh an d in the theory of Bose-Einstein condensation 
|12l[l^[l4|. However, recently it has been suggested that 
water-wave GSs are possible too, over periodic bottom 
boundary pj| In the cited works [H, Ojj], planar 

potential flows with one-dimensional (ID) free bound- 
ary were studied both analytically — with the help of 
an approximate model possessing relatively simple par- 
ticular solutions, and numerically — exact equations of 
motion for an ideal fluid with a free surface over a non- 
uniform bed were simulated, in terms of the so called 
conformal variables (see Ref. [17]). For two horizontal di- 
mensions, the question about possibility of GSs or some 
other coherent water-wave structures over periodic bed 
profiles is still open. The present work is a step to study 
this problem. More specifically, we suggest here approx- 
imate equations of motion for a two-dimensional (2D) 
free surface over a quasi-one-dimensional, locally peri- 
odic non-uniform bed. Then we present some numerical 
results where spontaneous formation of localized nonlin- 
ear structures is clearly seen, and the process looks as 
a kind of wave collapse. It should be emphasized that 
the Bragg interaction between the forward wave and the 
backward wave plays the key role in this process, since the 
wave dynamics over a flat bottom is rather different, and 
the extreme waves there are not so high as they are with 
the Bragg interaction. It is quite possible that effects 
considered in the present work have in some cases rela- 
tion to the extensively discussed phenomenon of rogue 
waves i Eil ! namely at seas with a non-uniform depth 
~ 50 - 100 m. 

Analytical and numerical results of Ref. [161 ] imply that 
the most suitable asymptotic regime for observation of 
nonlinear Bragg structures of the GS type in water- 



wave systems is the relatively short-wave (or deep-water) 
regime, when an effective water depth ho and the main 
Bragg-resonant wave number k form a small parameter 
e = cxp(— 2nho) ~ 0.01. It should be noted here that 
the opposite long-wave asymptotic regime is very differ- 
ent, and typical nonlinear structures at shallow water are 
KdV-type soiitons (see, Refs.[20l l2l|). In the considered 
deep-water regime, one can use the fourth-order Hamil- 
tonian functional 7i{?7, ip} for weakly-nonline ar g ravity 
water waves in the approximate form (see Ref. [lq] ) . 
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where r = (x, y) is the position in the horizontal plane, g 
is the gravity acceleration, the canonical variables 77(1", t) 
and ip(r,t) are the vertical coordinate of the free sur- 
face and the boundary value of the velocity potential 
respectively, the linear operator k = (k 2 + k 2 ) 1 / 2 is di- 
agonal in the Fourier representation, and K = k + B 
is a linear operator connecting the velocity potential 
ip(x,y,z = 0) at the unperturbed free surface and the 
quantity d z ip(x,y, z = 0), in the presence of an inhomo- 
geneous bed. The operator B depends on a given bed 
profile in a complicated manner (see Refs.pj| [22j], where 
some expansions of this operator are discussed), but it 
is definitely "small" in the deep-water limit, B ~ e, and 
that is why we keep B only in the quadratic part of the 
Hamiltonian, while in all higher-order terms we write k 
instead of K. 

Canonical equations of motion, corresponding to 
Hamiltonian |T]), are still too complicated to be treated 
analytically, and very time-demanding when being solved 
numerically. Therefore we shall consider here a simplified 
model which can be derived from |T]) in some limit. To 
do this, let us recall the well known fact that a suitable 
weakly nonlinear canonical transformation, 
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(where lu^ = (glkQ 1 / 2 is a linear dispersion relation in the 
absence of the bottom boundary) can exclude the third- 
order terms from the deep-water Hamiltonian, as well as 
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the non-resonant wave interactions <-> 4 and 1 <-> 3 (see 
Refs.JH, Hil])- Accordingly, the shape of the free surface 
r)(x,y,t) is given by the formula below: 
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In terms of the new normal complex variables 6k (t), the 
wave dynamics in deep-water regime is described by the 
following integral equation, 

ib k as uj k b k + Lb k + . . . + i J T(k, k 2 ; k 3 , k 4 )b k2 b k3 b k4 
x<5(k + k 2 - k 3 - k 4 )dk 2 dk 3 dk 4 , (4) 

where L is a "small" linear non-diagonal operator related 
to B, and T(ki, k 2 ; k3, k 4 ) is a known continuous func- 
tion (but the corresponding explicit expression is rather 
complicated, see Refs. [23l. |24|). For our purposes it is 
sufficient to know that in purely one-dimensional case 

T(k 1 ,k 2 ;k 1 ,k 2 ) cx |fci||fc 2 |(|fci +k 2 \ - \h-k 2 \), (5) 

and thus T(k, -k; k, -k) = -T(k, k; k, k) < 0. 

Let us assume for the moment that a bed profile is 
strictly x-periodic, with the period A = tt/k, and a ID 
wave spectrum is concentrated near the Bragg-resonant 
wave vectors ±k = ±ne x . Now we introduce two 
slowly- x-dependent functions, the envelopes A±(x,t) of 
the forward- and backward-propagating waves, 

rj(x, t) = Re[A + exp(inx — iujt) + A_ exp(— inx — iu)t)], 

where uj = (gn) 1 / 2 (the corresponding deep-water wave 
period is T = In jib). Then we take into account only 
the main-order components of the integral operator L: 
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where e is a real positive number, which should be 
identified as e = cxp(— 2kHq), for some efficient depth 
ho, since the finite-depth dispersion relation is uj = 
[<7Ktanh(/io^)] 1 ^ 2 ~ &[1 — exp(— 2nho)]', a complex pa- 
rameter A = |A| exp(i(/>o) bears information about the 
width of the main Bragg gap (namely, the edges of the 
frequency gap are at u>i f2 w Q(l — e =F |A|), and about 
the phase 4>o of the main Fourier harmonics of the bed 
undulation in some suitable representation (see Ref.flU] 
for details). It is worth mentioning that for all ID bed 
profiles considered in the previous studies [HI, , the in- 
equality |A| < £ holds (|A| approaches e from the below 
if the bottom boundary consists of very narrow barriers). 

As a result of the standard procedure, we obtain the 
system of two coupled approximate equations for the 
wave envelopes A±(x,t) (compare to Ref.[16(), 
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where the dots mean the second- and higher-order x- 
derivatives. The above equations describe weakly nonlin- 
ear ID water-wave GSs quite well: known solitary-wave 
solutions of this system were compared in Ref. [la ] to nu- 
merical simulations of the exact equations of motion for 
an ideal fluid with a free surface, and a reasonable corre- 
spondence was found. 

Now it is clear how to generalize the above system to 
the 2D case: we have just to add ^-dependent dispersive 
terms to the left hand sides, which appear from expan- 
sions of the operators [— k~ 1 / 2 ((— id x ± k) 2 — d 2 ) 1 ^ + 1]. 
Moreover, we shall consider here the case when s and 
A are no longer constants, but they are some functions 
slowly depending on the horizontal coordinates x and y. 
This assumption makes our model more realistic and rich, 
since there is no perfectly periodic bed structures in the 
nature, while roughly periodic bars are quite often. Then 
the following natural generalization of Eqs. ([6][7j) can be 
suggested, 
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(^-^-h + J^ + -)A- = -e( X ,y)A. 
\ uj 2k ok 4kt / 

+A*(x, y)A + + ^ (|A_| 2 - 2\A + \ 2 ) A_, (9) 

where the dots mean omitted third- and higher-order par- 
tial derivatives. We keep here the second-order disper- 
sive terms because they are definitely important in the 
dynamics of relatively short wave groups. The linear and 
nonlinear terms in the right hand sides are of the same 
order of magnitude at ~ 0.1 and e ~ |A| ~ 0.01 

(let us note that water waves become strongly nonlinear 
if k|^4±| gl 0.3). The above equations have the standard 
Hamiltonian structure, iuj~ 1 dtA± = 5H,/5A* ± , with 



H = J[A* + D + A + + A*_D-A-]d 2 r 

+ J [AA* + A_ + A*A*_A+ - e(A* + A+ + A*_A_)]d 2 r 
+ k 2 ( [(|yl + | 4 + |A_| 4 )/4 — |A + | 2 |yl_| 2 ]<i 2 r, (10) 



(6) 



where the dispersive operators D± are 

D ± = T*0 b /(2k) + d 2 x /(8K 2 ) - d 2 J{AK 2 ) + ■■■ (11) 

It should be emphasized that a (spatially non-uniform) 
linear Bragg coupling between the forward wave and the 
backward wave introduces essentially new effects in the 
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FIG. 1: (Color online). The absolute value of A(x,y). 



dynamics, in comparison with the flat-bottom case. It 
is already seen from the linear- wave dispersion relations 
corresponding to the simplest case e = const and A = 
const > 0. There are two branches in the spectrum (type- 
1 and type-2 waves), behaving quite non-trivially, 
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Various nonlinear processes occur in this system, which 
create instabilities near corresponding resonant curves in 
k-plane. For example, two type-1 waves with ki = 
decay into two type-2 waves with wave vectors ±k 2 (2 <-> 
2 process), if 



2fii(0) « fi 2 (k 2 ) +n a (-ka). 



(13) 



Thus, a spatially homogeneous solution A + = A_ = 
Aq exp(— iflot), with a small but finite constant ampli- 
tude Aq, is unstable. A standard linear analysis shows 
that a maximum of the instability increment is reached at 
the perpendicular direction, near k 2 = ±e y ny/8A. More- 
over, besides the above indicated instability, there is an- 
other, long-scale modulation instability of type-1 waves, 
which completely surrounds the wave vector k = 0, un- 
like the situation for oppositely propagating waves at in- 
finitely deep water (see, e.g., Refs.[2a.l2q]). Therefore the 
tendency towards spontaneous formation of big waves is 
more strong if the bottom is (locally) periodic. However, 
a more detailed analytical study of Eqs.(0[9]), including 
full stability analysis and particular solutions, will be a 
subject of future research. In this work we only present 
results of numerical simulations which give a general im- 
pression about the dynamics, with particular attention 
to the phenomenon of wave collapses. 

The numerical simulations were performed in the stan- 
dard dimensionless square domain 2ir x 2tt, with periodic 
boundary conditions and with k = 100. In this example, 
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FIG. 2: The global maximum of wave amplitude versus time. 
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FIG. 3: (Color online). The absolute values \A±(x,y,t)\ at 
t/To ~ 150. The localized energetic structure is seen near the 
right upper corner of the \ A- \ map. 
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function e was taken in the form e = 0.014(1 + 0.1 cosx + 
0.1 cosy), while A(x,y) = 0.007[1 + 0.05£(a;, y)] was 
used, where a complex function £(x, y) contained discrete 
Fourier harmonics within — 5 < k x < 5 and — 5 < k y < 5, 
with amplitudes decaying as exp[— 0.1(£^ + ky)] and with 
pseudo-random phases (see FigfTJ where a map of the 
absolute value \A(x, y)\ is shown). The initial conditions 
were kA + = 0.08 + a + (x,y) and kA- — a_(x,y), where 
functions a±(x,y) had the Fourier spectra with ampli- 
tudes behaving as 0.001 exp[-0.04(fc^ + ky)] and with 
pseudo-random phases. Thus, the most part of wave 
energy was initially concentrated in the forward wave. 
However, after few tens of wave periods, the energy has 
been redistributed between the forward and backward 
waves due to the inhomogeneous Bragg coupling. Wide 
regions of slightly higher wave amplitude were formed, 
as a result of interplay between random initial conditions 
and random A, which process is reflected in Figf21 After 
that, self-focusing nonlinear interactions came into play 
and rapidly produced the big wave which is seen in FigJ3] 
It is interesting to note that the localized structure was 
developed only in one of the two wave envelopes: there is 
no clear sign of the same big wave in the opposite-wave 
amplitude. 

The x-size of the developed high wave group contains 
just 1-2 wave lengths, while the y-size is larger, about 10 
wave lengths. Strictly speaking, the assumption of nar- 
row spectra for A± is violated in this situation, so the 
system (|8]|9|) can describe the final stage of wave collapse 
only qualitatively. To get indirect confirmations for the 
wave collapses, we compared the above approximate re- 
sults to analogous simulations of full systems of the type 
(|4|), for some simple functions T m (ki, k2; k3, k4) which 
possess the same property (J5J) and thus mimic the true 
complicated matrix element T(ki, k2; k3, k 4 ). In partic- 
ular, we considered 

T m (X |ki|*|k 2 |*|k3|'|k4|"(|ki +k 2 | + |k 3 +k 4 | 



-|ki - k 3 | - |ki - k 4 | - |k 2 - ka| - |k 2 - k 4 |), (14) 

and there the wave collapses were observed as well (we do 
not discuss here subtle points related to such systems) . In 
the reality, however, higher-order nonlinearities become 
important at |k^4±| 0.3 and they produce wave break- 
ing manifested as the well-known "white caps". There- 
fore the question about maximal wave height at the finite 
stage of collapse can be fully answered only through real- 
world experiments. 

To conclude, in the present work we have derived and 
then numerically simulated nonlinear equations of mo- 
tion for complex water-wave amplitudes coupled both 
linearly through the (spatially non-uniform) Bragg in- 
teraction, and non-linearly through the cross-modulation 
terms. In the numerical experiments, we have observed 
spontaneous formation of highly energetic localized struc- 
tures which look like wave collapses rather than like soli- 
tons. The final stage of these wave collapses is actu- 
ally beyond applicability of our weakly-nonlinear model. 
That is why some real-world experiments or at least fully 
nonlinear simulations are desirable. Unfortunately, large- 
scale simulations in the framework of the exact Eule- 
rian dynamics, for instance with the help of the exist- 
ing boundary integral methods, are hardly possible at 
the moment because of their extremely time-demanding 
implementations. Perhaps, the weakly non-planar but 
fully nonlinear equations of motion for potential water 
waves over quasi-one-dimensional topography derived in 
Ref. [13] , might be useful as an intermediate step towards 
accurate numerical results. 
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